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PREFACE 


The  work  described  in  this  report  was  part  of  an  ILIR  project  (work  unit 
4A181101A91D01)  conducted  in  the  Center  for  Artificial  Intelligence  of  the  Research 
Institute,  USAETL.  Its  purpose  was  to  study  the  general  properties  of  fractals  and  to 
determine  some  of  their  potential  applications  to  mapping  and  the  topographic 
sciences.  Particular  goals  included  the  development  of  rudimentary  computer 
programs  for  generating  synthetic  terrain  and  the  elucidation  of  “fractal  dimension” 
in  terms  of  surface  roughness  and  natural  textures. 

During  the  period  of  the  investigation,  COL  Alan  L.  Laubscher,  EN.was  the 
Commander  and  Director,  Mr.  Walter  E.  Boge  was  the  Technical  Director  and  Dr. 
Robert  D.  Leighty  was  the  Director  of  the  Research  Institute. 

The  work  was  performed  in  the  Center  for  Artificial  Intelligence  and  the  report 
was  prepared  under  the  supervision  of  Lawrence  A.  Gambino,  Director,  Research 
Institute. 
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COMPUTER  GENERATION  OP  FRACTAL  TERRAINS 


INTRODUCTION 

The  concept  and  the  word  “FRACTAL”  were  introduced  by  Benoit  B. 
Mandelbrot  around  1975.  He  made  the  important  observation  that  many  objects  in 
the  natural  world  are  not  continuous,  smooth  or  measurable  by  the  traditional 
standards  of  the  applied  mathematician  but  require  the  use  of  more  general  methods 
and  constructs  of  pure  mathematics.  Furthermore,  they  are  characterized  by  a 
quantity  which,  for  any  particular  object,  remains  invariant  regardless  of  the  scale  of 
observation.  In  the  case  of  natural  objects,  this  quantity,  the  ‘fractal  dimension’, 
may  retain  its  invariance  over  many  orders  of  magnitude.  The  importance  of  the 
concept  can  be  surmised  by  the  fact  that  most  physical  theories,  aside  from 
fundamental  laws,  break  down  or  lose  their  validity  over  a  much  smaller  range  due  to 
inapplicability  of  the  underlying  assumptions. 

Since  their  introduction,  there  has  been  a  nearly  explosive  interest  in  fractals  in 
physics,  in  the  natural  sciences,  in  computer  graphics,  and  many  other  areas;  as  well 
as  for  their  existence  as  objects  that  are  beautiful  and  interesting  in  their  own  right. 
Several  of  the  most  well-known  ‘fractals'  have  a  long  existence  in  classical 
mathematics  and  have  been  historically  important  in  the  development  of  the  theory 
of  point  sets,  dimension  theory,  measure  theory  and  topology.  Continued 
developments  in  applications  of  fractals  will  probably  result  in  the  relocation  of  the 
superficial  borderline  between  applied  and  pure  mathematics  deep  into  the  territory  of 
the  latter. 

It  was  recognized  by  Mandelbrot  that  such  terrestrial  features  as  coastlines  and 
terrain  were  fractal.  As  such,  one  of  the  parameters  used  to  describe  a  given  instance 
is  a  ‘dimension’  that  need  not  be  (usually  is  not)  a  whole  number.  In  the  work 
described  here,  a  computer  algorithm  (due  to  Fournier,  Fussell  and  Carpenter)  was 
used  to  generate  fictitious  computer  representations  of  terrains.  The  display  of  the 
resulting  data  was  in  the  form  of  profiles  of  the  synthetic  terrain  surfaces.  The 
terrain  figures  to  be  presented  were  made  by  directly  photographing  the  video  display 
of  the  computer  terminal.  No  consideration  was  given  to  the  generation  of 
perspective  views  or  to  the  problems  of  representing  (“colored”)  surface  reflectances 
under  conditions  simulating  natural  lighting.  A  Gaussian  distribution  was  used  to 
generate  a  series  of  surfaces  for  which  the  fractal  dimension  varied  from  1.05  to  1.80. 
It  was  of  interest  to  determine  the  parameters  that  would  lead  to  realistic  simulations 
of  real  terrains. 

Another  important  application  of  fractals  to  the  topographic  sciences  is  the 
determination  of  the  fractal  dimension  for  given  real  natural  objects  and  the  use  of 
the  dimension  as  a  descriptor  or  identifier.  This  problem  is  called  the  inverse  problem 
and  is  not  treated  in  the  present  report. 


BASIC  FRACTAL  CONCEPTS 


Fractals  have  been  defined  formally  in  the  following  way. 

Definition.  A  fractal  is  a  set  for  which  the  Hausdorff  Besicovitch 
dimension  strictly  exceeds  the  topological  dimension. 

It  is  apparent  that  this  definition  contains  two  types  of  dimension. 

The  topological  dimension  PT  for  a  set  will  always  be  a  whole  number  and 
corresponds  with  our  intuitive  concept  of  dimension.  It  is  defined  by  assigning  the 
dimension  -1  to  the  empty  set  and  using  the  definition  below. 

Definition.  A  space  A'  has  dimension  <  Dr  if  every  point  p  can  be 
separated  by  a  closed  set  of  dimension  <  Dr  -  \  from  any  closed  set  not 
containing  P. 

In  this  way,  a  set  of  isolated  or  disconnected  points  has  DT  =  0  since  the  points  are 
already  disconnected  by  the  null  set.  Similarly,  ordinary  (non-fractal)  curves  are 
found  to  have  DT  =  l  since  any  point  can  be  separated  from  the  rest  of  the  curve  by 
removing  points.  And,  for  surfaces,  DT  =  2,  the  segmentation  is  accomplished  by 
removal  of  the  points  on  a  curve  Dr  =  l.  When  dealing  with  fractal  sets,  the  same 
principle  holds  but  much  greater  care  must  be  taken  because  of  the  complexity  of 
fractal  sets  and  sometimes  ingenious  constructions  may  be  required. 


On  the  other  hand,  the  Hausdorff  Besicovitch  dimension  D  is  a  measure  theoretic 
concept  and.  depending  upon  the  set,  its  determination  can  be  extremely  difficult. 
The  set  of  points,  P,  to  be  measured  is  covered  with  (closed)  balls,  ,  of  maximum 
radius  f  such  that 


n(<) 

r  c  (j/?,(<) 

1*0 

The  measure  of  a  k-dimensional  ball  of  radius  /?  is  given  by 


Mk(m  =  «(*)  Hk 


2 


and  where  a  is  given  in  terms  of  the  f  function. 


o(t)  =  - - 

The  Hausdorff  measure,  of  dimension  A.  of  the  set  I'  is  then  given  by 

Hk(P)  =  lint  inf  U(A)  ,r.(Oi* 

—o*  l  ,_o 


and  where  t  he  lim  inf  is  to  he  taken  over  all  possible  coverings.  For  ordinary  (non¬ 
fractal)  sets,  k  is  taken  as  the  dimension  0T  and  the  Hausdorff  measure  for  any 
reasonable  case  turns  out  to  be  the  corresponding  count,  length,  area,  or  volume.  For 
fractal  sets,  k  need  not  be  an  integer  and  there  will  he  one  particular  value.  k=D,  for 
which  the  Hausdorff  measure  is  positive  and  finite. 

()•  Hd  *  oc 

This  dimension.  D.  is  called  the  Hausdorff  Besicovitch  dimension  and  is  the  most 
fundamental  parameter  defining  the  fractal.  The  measure  itself  is  usually  of  little 
importance  since  it  is  difficult  to  give  it  a  physical  meaning  for  non-integral  D.  For 
k>D  we  will  have  //*(P)  =  0  and  for  k  <  D  the  corresponding  measure  is  infinite 
//‘(P)  =  oc.  It  is  well  known  [Hurewicz  &  Wallman]  that  the  Hausdorff  Besicovitch 
dimension  is  never  less  than  the  topological  dimension.  However,  there  are  fractals 
with  dimension  given  by  an  integer. 

Because  of  the  complexity  of  the  above  concepts,  a  more  appropriate  working 
definition  has  been  adopted  for  the  present  study. 

Definition.  Fractals  are  mathematical  constructs  that  give  realistic 
representations  of  natural  patterns  and  textures  and  that  exhibit 
invariant  properties  under  scale  change  (magnification). 

For  geometrically  regular  fractals,  enlargement  by  a  certain  factor  will  bring  a  subset 
of  the  fractal  set  back  into  coincidence  and  this  accounts  for  the  striking  appearance 
of  many  fractals.  For  the  surfaces  to  be  generated  here,  the  fractal  should  exhibit  the 
same  statistical  properties  under  enlargement.  This  is  insured  by  the  generating 
algorit  Inn. 
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Equally  important  in  this  study  is  the  question  of  whether  the  fractal  has  an 
appearance  similar  to  that  of  a  natural  terrain.  Specifically,  what  dimensions  and 
other  generating  parameters  lead  to  realistic  results. 


THE  FRACTAL  GENERATION  ALGORITHM 


The  method  to  be  described  for  generation  of  fractal  surfaces  was  originated  by 
Fournier,  Fussell  and  Carpenter1  and  is  usually  known  as  the  FFC  algorithm, 
although  this  designation  may  confuse  it  with  a  different  algorithm  also  given  by 
them.  This  algorithm  was  chosen  because  of  the  simplicity  of  the  computation,  its 
flexibility  and  adaptability  to  potential  applications  in  the  topographic  community, 
its  realism,  and  its  general  compatibility  with  “standard”  methods  of  computer 
graphics.  These  criteria  are  discussed  in  the  next  paragraph. 

The  FFC  algorithm  to  be  described  allows  for  local  computation  of  terrain 
elevations  by  densification  of  points  from  a  set  of  elevations  given  over  a  square  grid. 
The  computation  requires  only  the  averaging  of  elevations  at  the  corners  of  a  square 
and  the  addition  of  a  scaled  random  number  (according  to  a  fractal  generating  rule). 
This  permits  filling  of  fictitious  terrain  textures  in  cases  where  actual  topographic 
data  may  be  given  on  a  sparse  (coarse)  grid.  It  also  permits  the  filling  of  more  detail 
in  regions  of  higher  magnification,  such  as  in  the  foreground,  or  in  .cases  of 
“zooming-in."  The  method  also  allows  for  variation  of  fractal  dimension,  vertical 
dilation,  or  the  distribution  function  from  one  local  area  to  another.  Over  square 
regions  where  the  grid  size  does  not  vary,  only  elevations  need  be  stored  since 
horizontal  coordinates  are  easily  generated  from  grid  indices  and  from  the  grid 
spacing  and  origin. 

Consider  the  situation  shown  in  figure  la.  The  grid  lies  in  a  horizontal  plane 
and  elevations  are  known  for  the  points  given  by  the  intersections  of  the  grid  lines 
(including  those  on  the  bounding  sides).  Elevations  are  to  be  generated  over  the 
center  point  of  each  square  as  depicted  here  by  the  dots. 


i 

FIGURE  1.  FRACTAL  FILLING  OF  THE  TERRAIN  GRID. 

1  Alain  Fournier,  Don  Fussell.  and  Loren  Carpenter.  "Computer  Rendering  of  Stochastic  Models"  (  omiminir.nlloii  of 
the  ACM,  vol  25,  no  6,  June  1982,  pp.  371-384 


For  each  small  square  cell,  the  average  elevation  E„  of  the  corner  (vertex)  points 
is  taken.  A  pseudo-random  number  chosen  from  a  Gaussian  distribution  R  is  used  to 
generate  the  new  elevation,  £«„,  the  center  point  by 


=£..+/?  P-D 


where  6  is  the  distance  from  the  center  point  to  the  vertex.  Here  D  is  the  fractal 
dimension  and  3  is  the  dimensionality  of  the  embedding  space. 

After  this  has  been  repeated  for  each  cell,  the  elevations  will  be  known  over  a 
finer  diagonal  grid  containing  both  the  original  elevations  and  the  newly  generated 
ones  as  shown  in  figure  lb.  The  procedure  is  repeated  again  for  elevations  of  the 
center  points  of  this  grid  with  the  exception  that  for  outer  bounding  edges,  the 
average  cell  elevation  is  determined  from  only  three  points. 

At  this  stage  a  new  (horizontal-vertical)  grid  has  been  established  as  shown  in 
figure  lc.  It  will  contain  four  times  as  many  points  as  the  initial  grid  in  figure  la. 
The  complete  procedure  can  then  be  repeated  ag'Mn  and  again  until  a  sufficiently 
dense  set  of  elevations  has  been  computed. 

After  the  elevations  have  been  computed,  a  vertical  dilation  is  applied  to  adjust 
the  amount  of  vertical  exaggeration.  It  was  found  that  surfaces  of  larger  fractal 
dimension  required  much  less  shrinkage  in  the  vertical  direction  than  those  with 
fractal  dimension  near  2. 
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EXAMPLES  OF  FRACTAL  TERRAINS 


Examples  of  computer  synthesized  terrain  are  given  in  appendix  A.  They 
include  results  for  various  dimension  and  vertical  exaggeration  and  should  serve  to 
demonstrate  the  effects  of  varying  these  parameters.  The  distribution  that  was  used 
in  all  cases  was  an  approximated  cumulative  distribution  function  for  a  Gaussian. 
The  views  of  the  surface  displayed  on  the  screen  are  not  in  perspective,  but  are  drawn 
in  the  same  orthographic  way  that  three-dimensional  mathematical  objects  are  drawn 
(e.g.  on  a  blackboard).  The  three-dimensional  appearance  is  partly  an  illusion  since 
only  profiles  of  the  surface  are  plotted.  Nevertheless,  since  most  people  are 
accustomed  to  such  plots,  they  can  give  an  excellent  feel  for  the  roughness,  texture, 
and  other  qualities  of  the  surfaces  involved. 

The  effect  of  varying  the  amount  of  vertical  exaggeration  is  shown  in  terrains 
1,2,  and  3.  These  three  examples  all  have  fractal  dimension  D=  2.05  and  were 
produced  by  the  same  set  of  random  numbers.  They  were  produced  using  vertical 
dilation  factors  of  0.05,  0.1,  and  0.2,  respectively. 

The  computer  renderings  in  terrains  4,  5,  6,  7  and  8  show  the  effect  of  varying 
the  fractal  dimension.  The  dilation  factor  was  also  different  for  each  example.  For 
these  surfaces,  the  fractal  dimension  D  took  the  values  2.05,  2.1,  2.2,  2.4  and  2.8, 
respectively,  while  the  amount  of  vertical  exaggeration  for  each  case,  in  turn,  was 
0.05,  0.1,  0.2,  0.4  and  0.8.  The  difference  in  these  two  parameters  is  clearly  seen  by 
comparing  terrains  4  and  8.  The  former  is  much  smoother  than  the  latter,  which 
shows  roughness  much  greater  than  any  imaginable  real  world  terrain,  while  the  total 
vertical  exaggeration  is  not  drastically  different.  Terrains  5  and  6  are  more  rugged 
than  4  and  show  greater  vertical  differences  than  either  4  or  8,  while  remaining  much 
smoother  than  terrain  8. 

An  example  of  an  undesirable  result  is  shown  in  terrain  9.  The  unnatural  spikes 
on  the  back  edge  are  believed  to  result  from  the  chance  occurrence  of  choosing 
numbers  far  out  on  the  tails  of  the  generating  Gaussian.  This  can  be  eliminated  by 
using  a  different  distribution  that  does  not  possess  such  long  tails;  for  example,  a 
chopped  off  Gaussian.  The  rarity  of  such  occurrences  may  be  surmised  by  their 
absence  from  the  other  renderings  since  each  one  represents  263,169  generated 
elevations  on  a  513  x  513  grid  with  each  8th  profile  plotted. 


7 


THE  LISP  COMPUTER  PROGRAMS 


A  listing  of  the  computer  programs  is  given  in  Appendix  B.  They  were  written 
in  the  LISP  language  for  compatibility  with  other  work  in  the  Center  for  Artificial 
Intelligence,  ETL,  and  because  of  the  available  hardware.  A  brief  description  of  the 
principal  functions  will  be  given  here. 


Before  running  a  program  to  generate  a  fractal  surface,  the  function  foooz  must 
be  executed  to  initialize  the  program  and  to  define  some  of  the  arrays  needed.  The 
array  nordis  contains  a  numerical  table  representing  the  cumulative  distribution 
function  for  a  Gaussian.  The  function  u-to-g  converts  a  random  number  from  a 
uniform  distribution  to  a  random  number  chosen  from  a  Gaussian  distribution. 

The  fractal  generated  is  three  dimensional,  but  the  indices  of  a  two-  dimensional 
array  are  used  for  horizontal  coordinates  while  the  stored  numbers  give  the  vertical 
coordinate.  The  functions  ssu  and  ssv  generate  screen  coordinates  from  given  t  hree- 
dimensional  coordinates.  They  do  not  give  perspective  views  but  only  a  convenient 
readily  visualized  representation.  They  are  used  by  dsl  to  draw  a  straight  line  given 
its  three-  dimensional  endpoints. 

The  function  elev  generates  a  fractal  elevation  in  the  center  of  a  cell  in  the 
fractal  grid.  Its  arguments  are:  the  average  elevation  of  the  corner  points,  the  fractal 
dimension,  and  the  distance  from  the  center  of  the  cell  to  a  corner. 

The  function  surfract  generates  the  fractal  surface  by  using  the  previously 
described  filling  algorithm.  Its  sole  argument  is  the  fractal  dimension.  Some 
intermediate  views  are  displayed  on  the  screen  to  monitor  the  filling  operation. 

The  function  dilate  is  used  to  adjust  the  amount  of  vertical  exaggeration  by 
applying  a  dilatation  factor.  The  arguments  of  this  function  are  the  scaling  factor 
and  another  number  specifying  how  many  profiles  should  be  skipped  before  plotting 
the  next  one.  It  is  interesting  to  note  that  a  negative  scaling  factor  has  the  effect  of 
inverting  the  surface,  i.e.  exchanging  mountain  peaks  and  valleys. 

The  function  stratum  is  used  to  put  a  base  on  the  surface  and  to  give  it  the 
appearance  of  a  solid  slab.  Its  first  two  arguments  specify  the  thickness  of  the  slab 
and  the  distance  between  profiles,  respectively. 

Typical  running  times  for  a  513  x  513  surface  generation  and  display  have  been 
a  little  less  than  10  minutes. 


I 


CONCLUSION 


Computer  programs  were  written  in  LISP  that  make  use  of  a  Fournier,  Fussell, 
Carpenter  algorithm  to  generate  synthetic  fractal  terrains  having  various  dimensions. 
It  was  found  that  the  degree  of  realism  was  best  for  fractal  surfaces  with  dimension 
closer  to  2  and  that  those  with  dimension  closer  to  3  were  much  too  rough  for  any 
imaginable  real  world  terrains. 

It  is  believed  that  this  algorithm  can  be  used  for  filling  scenes  where  actual 
elevations  are  given  on  a  sparse  grid  in  order  to  achieve  realistic  rendering  of  the 
terrain.  To  accomplish  this  will  require  the  generation  of  perspective  views  where 
color  or  gray  levels  are  used  to  depict  surface  reflectances  under  simulated  real  world 
illumination  conditions.  Such  renderings  would  involve  considerably  more  research 
and  programming. 
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THE  LISP  COMPUTER  PROGRAMS 


21 


SM  6  »H«  111*  (rc.lip  P«««  I 


(4«lM  foe  at  (> 

(Hi  *b*5*  10.) 

(**t  *ibite  10.) 

toot  'wrl  (atka-array  *($U  513)  *:tyi>*  *ert-i  ' ;  initial-va  luo  0.0)) 

(•of  'saurt  (eeke-array  '(51)  513)  * t * yr  »  *ert-q  * : Ini tlal-walue  0.0)) 

(•ot  'jralm*  (tikfirrip  *(40  40)  *:ty|:e  *er  t-o  ’ :  initl  al-«aluo  0.0))) 

(aot  *nordn  (take -array  *(*1  2)  'Jtyna  *arl-e)) 

(f  Ularrey  nerdis  *(0.00  0.500  0.05  0.72}  9.10  O.St'D  0.15  0.5*0  0.20  0.5)0 
0.25  0 . >0 1  9 .  )0  0.410  0.3S  0.43)  0.40  0.451  0.45  0.4)4  0.50  0.401  0.55 

0.)-)0  0.40  0.  )24  0.45  0.142  0.20  O.fSt  0.)5  0.))3  0.00  0.  )00  0.05  0.002 

0.90  3.316  0.05  0.020  1.00  0.041  1.35  0.H5J  1.10  0.044  1.15  0.0)5  1.20 

0.005  1.25  0.304  1.202  0.000  1.30  0.943  1.35  4.011  1.40  0.010  1.45 

0.024  1.50  0.03)  1.55  0.9)9  1.40  0.I4S  1.445  0.050  1.49  0.051  !.)0 

0.055  1 . ) 5  0.040  1.00  0.0*4  1.05  0.3*0  1.30  0.0)1  1.95  0.014  1.140 

0.9)5  2.00  0.0)1  2.05  0.190  2.10  0.302  2.15  0.004  2.2*  0.000  2.25  0.090 

2.  JO  0.  )4*l  2.324  0.990  2.35  0.991  2.40  0.052  2.45  0.993  2.50  0.194 

2.524  0.035  2.45  0.094  2.)5  0.00)  2.90  0.392  3.00  0.900  3.00  1.100)) 


(Oafuii  u-to-.jia) 

(cond  ((Imp  a  0.5) 

(4«t  *•  -1.0) 

(»et  *a  (-  1.0  a))) 

(I  (set  *•  *1.0))) 

(Uh  far  1  fra*  0  until 

(>  ( jr •(  nor  dll  1  1)  a) 
lo 

(«it  *J  1>) 

(cond  {(<  j  51) 

(mod  <(>•  <-  Caref  nardli  (*  J  I)  1)  a) 
(-  a  (aref  nordis  J  )>)) 

(*  i  (aref  nordu  j  0))) 

(t  (a  •  (aref  nerdli  (♦  }  1>  0)>)>) 
(l  (»  *  (aref  nerdli  J  0 ))))) 


(defun  iiu  («  y  z) 

(♦  53.  (»  (♦  1.5  a)  (•  .25  y)))) 


(defun  «i«  (■  y  i  tilt) 

(-  400.  (<  tilt  (♦  (•  .25  y)  (•  1.5  t>)>>) 


(defun  dll  (■*  ye  re  at>  yb  zb) 

(tend  tern  in  el -1  o  'Idraa-lint  (lla  (.no  »e  yr  .•’<)>  (fl«  (is*  »e  ye  ir  1.5)) 
(tie  (v,u  ab  yb  zb))  (M<  (is*  a‘»  *1'  zl>  1.5)))) 


(Oefun  e lev  (uh  d»a  rr) 

(♦  ewe  (s  (u-to-j  (//  (float  (rendin  IlC'M))  1-1000.)) 
(*  rr  <-  3.0  dia))))) 
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Sap  4  •*:**  1)05  frc.ltp  2 


(Mfwi  turfrnct  (dla) 

(tat  * turf  (ii»k<-»rrty  *(JU  SI3)  *!ty«|  '»r(-  |  *:lnlU«t*y<lM  0.0)) 

(Mt  'ad  go  512) 

(•*«  'dill  (✓/  adge  2)) 

(l««p  far  n  frta  t  until  (<  diit  1) 

4* 

(«t<  *rr  (//  (Hot  *d|t)  (fort  2.0))) 

(loop  for  1  froa  dltt  by  ado*  until  (>  I  512) 
do 

(loop  for  j  f pup  Plot  by  tdgo  until  (>  J  51*2) 
do 

(not  (ulav  (//  («  (ar af  yurt  (-  1  dltt)  t-  J  dl«t)) 

(•raf  turf  (*  1  dltt)  (•  j  dial)) 

(•rat  turf  (•  1  dltt)  (»  j  dltt)) 

(•rat  turf  (♦  1  dltt)  (♦  J  dltt)))  4.0) 

dip  rr) 
turf  1  j))) 

( laoo  for  t  from  0  by  dial  until  (>  1  $12) 
da 

C 1  top  for  J  from  (If  (addp  (//  1  dltt))  0  dltt)  by  tdga 
until  (>  j  512) 
do 

(cond  ((equal  1  0) 

(••at  (alt*  (//  (*  (nraf  tu'f  (>  1  dltt)  J) 

(araf  turf  1  (»  j  dltt)) 

(arat  turf  i  (-  J  dltt)))  ).0) 
dip  (float  dit«))  turf  1  j)) 

((•anal  1  512) 

(atot  (ala*  (//  <♦  (araf  turf  (•  1  dltt)  J) 

(araf  turf  I  (t  j  dltt)) 

(araf  turf  I  (-  j  dial)))  J.0) 
dip  (float  dltt))  turf  1  J)> 

( ( atual  J  0) 

(atat  (ala*  (//  (♦  (araf  turf  (t  1  dltt)  J) 

(araf  turf  ( -  1  dltt)  J) 

(araf  turf  I  (♦  J  dltt)))  J.0) 
dip  (float  ditt))  turf  1  J)) 

((oquat  J  512) 

(atat  (ala*  l//  (♦  (araf  turf  (♦  t  dltt)  J) 

(araf  turf  (-  1  dltt)  j) 

(araf  turf  I  (-  J  dltt)))  J.0) 
dip  (float  dltt))  turf  |  j)) 

(t 

(••at  (atey  (//  (♦  (araf  turf  (»  I  dltt)  J) 

(araf  turf  (-  I  ditt)  J) 

(tr»f  turf  1  (♦  J  dltt)) 

<  tr-f  t.,rf  l  (-  )  dltt)))  4.0) 
di a  (flo»t  dltt))  turf  1  J))))) 

(tat  'edge  ditt) 

(tat  'dltt  <2/  dltt  2)) 

(turf  44  turf)) 

(•and  tarninal-io  * : cl aar-alndoa ) 

(tarf  *4  turf)) 

(dafun  turf  (grata  sturf ) 

(leap  far  1  •ror  0  to  512  by  grata 
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|«»  *  0»»*»  III*  frc.lt*  III*  ) 


4* 

(loop  for  j  froa  1  to  >12 
do 

(dsl  (flppt  1)  C-  (fl«t(  j)  t.O)  («rtf  itpf  1  <-  j  t)> 
(float  1)  (float  J)  (trtl  *"vf  t  J)))) 

(U«l  fir  J  froa  0  Ip  >12  by  grata 
4p 

(loon  for  1  frpa  I  tp  >12 

Jo 

<4il  (-  (float  1>  |.q)  (float  j)  ( »r »f  iturf  (-  1  1)  J» 
(float  i)  (float  j>  ( or  if  taur f  1  J»>>>> 


(dafun  dllato  (scalp  dpi) 

(laa*  tor  1  ft  o  0  uni  1 1  >  1  >12) 

da 

(loop  for  J  fraa  0  until  (>  J  >12) 

ilo 

(asat  (•  scala  (arof  surf  1  J))  taurf  1  J))) 
(turf  dpi  saurf)) 


(do fun  sift  (scala  dol) 

(dllata  seal  >  dot) 

(laap  for  1  froa  0  until  (>  1  H) 
da 

(loop  for  J  fraa  •  until  (>  J  It) 
do 

<  »s#t  (arof  saurf  (•  6  (a  12  1))  (»  *  (•  11  j>)>  fratln«  1  j>>)> 


(dafun  stratua  (tMehnass  grata  (lopttanal  (saurf  saurf))) 

(dal  0.*  0.0  (-  0.0  tMehnass)  >12.0  0.0  (-  0.0  tMehnass)) 

(dsl  >12.0  (-.C  (-  0.0  tMehnsss)  >12.0  >12.0  (-  0.0  tMehnass)) 

(laap  for  1  froa  0  to  >12  by  grata 
do 

(dsl  (float  1)  0.0  (-  0.0  tMehnass)  (float  1)  0.0  (arof  saurf  1  0)) 

(dsl  >12.0  (float  1)  (-  0.0  tMehnass)  >12.0  (float  1)  (arof  saurf  >12  1)))) 
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